"Chiral" universality class behavior of a 
non-chiral antiferroquadrupole system 

M. Zukovic a *, T. Idogaki b , and K.Takeda a,b 

a Institute of Environmental Systems, Kyushu University, Fukuoka 812-8581, Japan 

b Department of Applied Quantum Physics, Kyushu University, Fukuoka 812-8581, Japan 

The planar Heisenberg system with antiferroquadrupolar exchange on 3D stacked triangular 
lattice is shown to belong to the new "chiral" universality class, predicted for chiral antifer- 
romagnets. The present system, however, displays no such chirality, which, according to the 
currently widely accepted concept, is supposed to be a key ingredient for the new critical 
behavior. Our claim that the new universality class should not be limited only to the chiral 
antiferromagnets is based on a simple mapping between the chiral antiferromagnetic and 
antiferroquadrupolar systems and supported by actual calculations of the transition tem- 
peratures and critical indices from finite-size scaling analysis of data issued from Histogram 
Monte Carlo simulations. In order to demonstrate the existence of two different universality 
classes in the behavior of quadrupolar systems, the finite-size scaling is also performed for 
the system with ferroquadrupolar exchange, which, in contrast to the antiferroquadrupolar 
exchange system case, produced standard critical behavior, as could also be anticipated from 
the mapping. 
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It has been theoretically predicted^ and supported by experiment^ that geometrically 
frustrated XY antiferromagnets on a stacked triangular lattice (STL), described by Hamil- 
tonian 



where Si = (Si x , Si y ) is a two-dimensional unit vector at the ith lattice site , denotes 
the sum over nearest neighbors, and J < is the bilinear exchange interaction constant, 
should belong to a new universality class linked to the two-fold chiral degeneracy inherent 
to the 120° ordered spin structure. Critical behavior of these systems is characterized by 
new critical indices, different from those for non-frustrated systems with XY spin symmetry. 

Motivated by the existence of two different universality classes in the critical behavior 
of the geometrically frustrated chiral antiferromagnets and non-frustrated systems with XY 
spin symmetry, we investigated critical behavior of the corresponding systems with only 
biquadratic exchange interaction. Namely, we studied models described by the Hamiltonian 



where J' is the biquadratic exchange interaction constant, in which we considered both 
J' > and J' < 0. For the case of J' > (ferroquadrupolar exchange), the system is 
non-frustrated and is expected to belong to the standard 3D XY universality class. On the 
other hand, the system with J' < (antiferroquadrupolar exchange) is frustrated due to the 
triangular lattice geometry, resulting in the non-collinear ground state (Fig.l). However, 
such a system, compared to the one with the antiferromagnetic exchange (Fig. 2), lacks the 
chirality - an essential element which is supposed to be responsible for the new critical 
behavior of the STL antiferromagnets. The question is: Will such a system still display the 
standard universality class behavior? 

Since the universality class is in principle determined by the space dimensionality and 
symmetry of the order parameter, let us first consider the issue of the order parameter 
symmetry in the cases of the STL antiferromagnet described by Hamiltonian ([1]) and the 
STL antiferroquadrupolar system described by Hamiltonian (J2]) (the space dimensionality 
is the same). XY continuous symmetry is common for both systems, however, there is 
an additional ingredient which makes the two systems different. In the case of the STL 




(1) 




(2) 
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antiferromagnet one may define a local chirality at each elementary triangular plaquette, by 



where the summation runs over the three directed bonds surrounding each plaquette, p, 
which is an Ising-like quantity representing the sign of rotation of the spins along the three 
sides of each plaquette. Due to the chirality the STL antiferromagnet has two-fold degener- 
acy of the ground state (k = +1 and k = — 1), resulting in the structure with spins arranged 
on plaquettes with turn angles +120° and —120°, respectively. A minimum energy condition 
is realized by an arrangement in which the + and — plaquettes alternate, as shown in Fig. 2, 
producing long-range chiral order at low temperatures. On the other hand, as far as the 
chirality is concerned, the STL antiferroquadrupolar system has four-fold degeneracy in the 
ground state of each plaquette (k = ±1,±|), resulting in the structure with four possible 
turn angles between two neighbouring spins ±120°, ±60° (Fig.l). However, unlike in the 
chiral case, there is no energetically favorable arrangement among the four kinds of plaque- 
ttes in the antiferroquadrupolar system and, hence, the plaquettes do not order even at low 
temperatures. So, compared to the chiral antiferromagnet with both Ising and rotational 
symmetries of the order parameter (Z 2 x 0(n=2)), the order parameter symmetry of the 
system with antiferroquadrupolar exchange is apparently different. In the following we will 
show, by a simple mapping between the dipolar and quadrupolar forces^, that this difference, 
however, should not constitute grounds for different critical behavior. 
The Hamiltonians H and %' can be rewritten in the respective forms 



where (pi — (pj is an angle between the neighbouring spins Si and Sj. Using the identity 
cos 2 (a) = (1 + cos2a)/2 we can write the corresponding partition functions in the forms 




(3) 




(4) 



and 




(5) 



Z 



( j; t ) = / n exp t cos< ^ ~ 





(6) 
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and 

/ 2J'\ N f 2w rr ( J' 

Z'{J',T) = (^Trexp — J J[dp k exp I — E cos2 (^ - (p d ) | , (7) 



It can easily be seen that if we substitute i?j — ■dj = 2(</?j — <pj) into Eqn. ([7]) and assume 
| J\ = | J' |, we have 

Z(J, T) = I 2tt exp - Z'(J', T/2) , (8) 

which holds for the case of a non-frustrated systems with collinear ordering, as well as the 
present frustrated systems case with non-collinear ordering. Therefore, in either case we 
can expect the same critical behavior of the systems ([T]) and ([2]) at the respective transition 
temperatures T c and T q = T c /2. 

In order to confirm the prediction based on the mapping, we further performed a His- 
togram Monte Carlo (HMC) simulation^ analysis. We considered the systems of N = L 3 
spins, where the lattices linear dimension L = 12, 18, 24 and 30, assuming periodic boundary 
condition throughout. Spin updating followed the Metropolis dynamics and averages were 
calculated using 2 x 10 6 Monte Carlo steps (MCS) after discarding another 1 x 10 6 MCS for 
thermalization. Then we used the data issued from the HMC simulation in order to perform 
a finite-size scaling analysis of the following physical quantities: the quadrupole long-range 
order (QLRO) parameter q, 



q = (Q) /N, where Q 



(9) 



the corresponding susceptibility per site Xg 



X « = Nk B T ' (10) 
and the logarithmic derivatives of (Q) and (Q 2 ) with respect to K = 1/ksT 

D '^M la{Q2) -^wr- {E) - < 12) 

Temperature-dependences of these quantities display extrema at the L-dependent transition 
temperatures T q (L), which at a second-order transition are known to scale with a lattice 
size as: 

Xq ,max{L) (X L^" , (13) 




L l/u q ^ 



(14) 



D 



(L) oc L 1 /"' 



(15) 



where u q and 7 g represent the correlation length and susceptibility critical indices, respec- 
tively. 

In Fig. 3 we present the scaling analysis for the case of J' > in ln-ln plot. As expected, 
the transition is of second order and the obtained values of the critical indices, v q = 0.667 
±0.008 and 7 g = 1.333 ± 0.018, are close to those for the standard XY universality class 
(see Table 1). Similar scaling analysis was done for the system with the antiferroquadrupo- 
lar couplings (J' < 0) and is presented in Fig. 4. The transition is again of second order, 
however, the values of the critical indices v q =0.520 ± 0.003 and 7 9 = 1.072 ± 0.009 turn 
out to be quite different from the standard ones but strikingly close to those of the new 
chiral universality class. We note that for the sake of comparison we also performed scaling 
analysis for the antiferromagnetic system ([1]) with J < and the calculated values of the 
critical indices, v = 0.514 ± 0.007 and 7 = 1.074 ± 0.013, (see Table 1 for comparison with 
the values obtained from some other studies) were indeed found to be in a nice agreement 
with those for the antiferroquadrupolar system. The comparison of the critical temperatures 
T c and T q also shows that they indeed differ by the factor two and, hence, along with the 
critical indices (all summarized in Table 2) excellently corroborate the prediction based on 
the mapping. Therefore, the present result would suggest that not necessarily the chirality 
as defined in equation fl3]) but rather frustration, resulting in the non-collinear ground state, 
is the element causing the novel critical behavior. 

In summary, the results presented in this work provided support for existence of the new 
universality class predicted by Kawamurar-^. However, the completely new result is that 
the system with antiferroquadrupolar exchange on STL, which has no chirality possessed 
by its counterpart with antiferromagnetic exchange, was found to display critical behavior 
strikingly similar to that of the chiral STL antiferromagnet of the new universality class. 
Based on the mapping between the two systems and HMC numerical results, we conclude 
that the two systems should indeed belong to the same universality class, that so far has 
been thought to be limited only to the chiral antiferromagnets. 
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TABLE I: Critical indices v, 7 for dipole ordering - previous studies. 





V 


7 


Ref. 


Standard XY model class 


0.669 


1.316 


10 


Chiral XY model class 


0.54 ± 0.02 


1.13 ± 0.05 


1. 


0.57 ± 0.03 


1.10 ± 0.05 


3 


0.54 ± 0.03 


1.01 ± 0.08 


4 
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TABLE II: Critical indices v q , j Q and v, 7 for quadrupole and dipole ordering, respectively. 





Vq 


7g 


T q /\J'\ 


J' > 


0.667 ± 0.008 


1.333 ± 0.018 


1.520 ± 0.002 


J' < ot 


0.520 ± 0.003 


1.072 ± 0.009 


0.729 ± 0.002 




V 


7 


T C /\J\ 


J < 


0.514 ± 0.007 


1.074 ± 0.013 


1.458 ± 0.002 



t Since the sign of the inter-plane interaction J' z is irrelevant in this case, for simplicity we used 

J' z >0. 
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FIG. 1: Ground state spin configuration of the system with antiferroquadrupolar exchange (J' < 0) 
shown in one triangular plane. In contrast to the frustrated XY antiferromagnets (Fig.2), the 
elementary triangles do not show chiral ordering. 
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FIG. 2: Ground state spin configuration of the system with antiferromagnetic exchange (J < 0). 
The "up" (+) and "down" (— ) triangles consisting of spins with turn angles 120° and —120°, 
respectively, show chiral ordering. 
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FIG. 3: Scaling of maxima of the quantities \q, D\ q and Di q in ln-ln plot for the case of J' > 0. 
The slopes yield values of \jv q for D\ q , D2 q and Jq/vq for \q- 
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